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ABSTRACT 

We consider stationary, forced, imbalanced, or cross-helical MHD Alfvenic turbulence where the 
waves traveling in one direction have higher amplitudes than the opposite waves. This paper is 
dedicated to so-called strong turbulence, which cannot be treated pcrturbatively. Our main result 
is that the anisotropy of the weak waves is stronger than the anisotropy of a strong waves. We 
propose that critical balance, which was originally conceived as a causality argument, has to be 
amended by what we call a propagation argument. This revised formulation of critical balance is 
able to handle the imbalanced case and reduces to old formulation in the balanced case. We also 
provide phenomenological model of energy cascading and discuss possibility of self-similar solutions 
in a realistic setup of driven turbulence. 

Subject headings: MHD - turbulence - ISM: kinematics and dynamics 



1. INTRODUCTION 

MHD turbulence appears in the dynamics of conduc- 
tive fluid in a generalized settings with large Reynolds 
numbers, or low physical dissipation. It is ubiquitous in 
the interstellar and intracluster medium, Earth magne- 
tosphere, solar wind, accretion disks, etc. In fact, it is 
laminar flows that constitute exception in astrophysics, 
while, generically, astrophysical fluids are turbulent. 

The study of MHD turbulence has been an old chal- 
lenge. First attempts to address it were classical papers 
by Iroshnikov (1963) and Kraichnan (1965) (henceforth 
IK model) . A good account for the state of the field could 
be found in Biskamp (2003). Usually turbulence is sub- 
divided into weak and strong, depending on the strength 
of non-linear interaction. While weak MHD turbulence 
allows analytical perturbative treatment (Ng & Bhat- 
tacharjee 1996, Galtier et al 2002, Chandran 2005), the 
progress in understanding strong turbulence came pri- 
marily from phenomenological and closure models which 
were tested by numerical simulations. 

Important theoretical works on strong MHD tur- 
bulence include Montgomery & Turner (1981), She- 
balin, Matthaeus, & Montgomery (1983), Higdon (1984). 
Those clarified the anisotropic nature of the energy cas- 
cade and paved the way for further advancement in the 
field. The study by Goldreich & Sridhar (1995, hence- 
forth GS95) identified the balance between perturbations 
parallel and perpendicular to the local direction of mag- 
netic field, i.e. "critical balance", as the key compo- 
nent of dynamics in strong magnetic turbulence. Al- 
though it dealt with incompressible MHD turbulence, 
GS95 also influenced further studies of compressible tur- 
bulence (e.g. Lithwick & Goldreich 2003). In particular, 
it identified the dominant role of Alfvenic perturbations 
for cascading of slow modes, which later confirmed with 
numerical simulations in both weakly and strongly com- 
pressive media (Cho & Lazarian 2002, 2003). 

Being a mean-field model in the spirit of Kolmogorov 
(1941), GS95 predicts the velocity and magnetic field 
fluctuation strengths and their anisotropies, in terms of 
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the local dissipation rate. Even though recently the sim- 
ple dynamical model of GS95 came under criticism (see, 
e.g., Boldyrev 2006, Gogoberidze 2007) with the motiva- 
tion to explain the deviations from GS95's -5/3 spectrum 
in numerical simulations (see, e.g., Muller, Biskamp & 
Grappin 2003), we feel that it does provide a good in- 
sight into MHD turbulence. The model that we present 
in this paper, is similarly to Kolmogorov or GS95, a mean 
field model, which does not account for any local dynam- 
ical effects or intermittencies (cf Beresnyak & Lazarian 
2006). These effects are beyond the scope of this paper 
and will be addressed elsewhere. 

While balanced MHD turbulence enjoyed much atten- 
tion, the opposite regime, i.e., imbalanced turbulence, 
was less developed. The analytical results were obtained 
for weak imbalanced turbulence (Galtier et al. 2002, 
Lithwick & Goldreich 2003) though they are applicable 
in a rather narrow range of imbalances. 

So far, the simulations of strong imbalanced turbulence 
were limited to rather idealized set ups (Maron & Goldre- 
ich 2001, Cho, Lazarian & Vishniac 2002), which did not 
allow making definitive conclusions about its properties. 
In this paper, in order to test our analytical model, we 
provide simulations which go a step further compared to 
the aforementioned studies. We hope that future higher 
resolution simulations will provide the definitive test for 
alternative models of imbalanced turbulence. 

We think that the best experimental data on the im- 
balanced regime is currently available from observations 
of solar wind turbulence (e.g., Horbury 1999). This data, 
collected by spacecrafts, is consistent with Kolmogorov 
—5/3 spectrum, but does not provide sufficient insight 
into the anisotropy with respect to the local magnetic 
field. The imbalanced turbulence is not a rare excep- 
tion, on the contrary, such processes as preferential decay 
of a weaker wave and the escape from the regions that 
generate perturbations make the imbalanced turbulence 
ubiquitous. It goes without saying that models of such 
turbulence are much needed in astrophysics. 

An attempt to construct the model of stationary 
strong imbalanced turbulence was done in Lithwick et 
al. (2007), henceforth LGS07. This model assumes the 
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strong GS95-type cascading of both large and small am- 
plitude oppositely moving modes. In what follows, we 
propose a different model of imbalanced turbulence. In 
§2 we revisit the critical balance argument and discuss a 
new process of cascading, which we relate to the process 
of propagating of Alfvenic perturbations in the field wan- 
dering which is induced by the oppositely moving wave. 
Using this "propagation cascading" in case of the bal- 
anced turbulence we recover GS95 relations. However, 
for the pronounced imbalance between the oppositely 
moving waves, this process results in a picture which is 
different from that in GS95. In §3 we discuss the scaling 
relations for the turbulent cascade that follow from our 
model and in §4 we compare our predictions with our 3D 
numerical simulations of the imbalanced cascade. In §5 
we compare our results with other studies. 

2. CRITICAL BALANCE REVISED 

The "perpendicular cascade" , a concept which was rig- 
orously developed in theory of weak Alfvenic turbulence 
(Ng & Bhattacharjee, 1996, Galtier et al, 2002), was a 
theory of nonlinear interacting Alfvenic waves which, due 
to the particular dispersion relation of the waves, con- 
served wave frequencies lo = k\\VA- It deeply contrasted 
with earlier Iroshnikov-Kraichnan models where, due to 
the assumed isotropization, the parallel wavenumber will 
be of the order of the total wavenumber and the fre- 
quency lo — k\\VA ~ kvA changes with k. The perpen- 
dicular cascade, however, makes the cascade stronger and 
not weaker, while the energy goes downscale. This raised 
a question of what happens when the perturbation theory 
breaks down and the turbulence becomes strong. GS95 
argued that the turbulence will stay on the edge of being 
strong, because of the uncertainty relation between the 
cascading timescale and the wave frequency t cas LO ~ 1. 
Their EDQNM closure model contained explicit ad hoc 
term that allowed for the increase of frequencies of the 
interacting wavepackets 1 . Since this process of increase 
of fc|| comes essentially from the irreversibility of energy 
cascading downscale, we call it the causality effect. This 
way nonlinear interaction stays marginally strong by con- 
trolling the anisotropy. But is this the only way to in- 
crease frequency or fc||? In this paper we advocate to 
supplement the causality effect with a different mech- 
anism which works when the cascading of the eddy is 
done by the counter-eddy, which is tilted with respect to 
the mean field due to the different definition of the local 
mean field at different scales. The details arc following. 

Indeed, the parallel wavenumber for an Alfvenic eddy 
in strong turbulence actually depends on how the local 
mean field is defined. Maron & Goldreich (2001) have 
shown that the Alfvenic eddy propagates along the field 
lines that are defined by the average field plus the field 
of the counter-wave. 

Let us try to calculate the characteristic uncertainty 
in fc|| that comes from this effect in the balanced case. 
In general, it will depend on eddy's polarization, but at 
maximum it will be around fcsin#, where 9 is the angle 

1 In this paper we will use wavevectors fcy, k± and length scales 
A = 1/feii, A = l/k± interchangeably. While wavevector repre- 
sentation highlights geometry and conservation laws, the length 
scales are measured quantities obtained from statistical averag- 
ing of numerical or observational data. We will also use "waves" , 
"wavepackets" and "eddies" interchangeably. 
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Fig. 1. — Upper: a w+ wavepackct, produced by cascading by 
w~ wavepacket is aligned with respect to w" wavepacket, but 
misaligned with respect to the local mean field on scale Ai, by the 
angle 8. Lower: the longitudinal scale A of the wavepackets, as 
a function of their transverse scale, A; A+, A - , Ai, A2 are the 
notations used in this paper. 



of the field wandering, or k±Sb(l)/vA (b is the magnetic 
field perturbation in the Alfvenic units). What is the 
characteristic scale I that enters this expression? We can 
argue that in the case of a local balanced turbulence there 
is only one designated scale, the one we consider cascad- 
ing from. Therefore, in the case of balanced turbulence 
we expect 8k\\ = k±Sb(l/k±)/vA- This coincides with 
the well known critical balance, derived from causality 
(indeed, in GS95 the cascading rate is of the order of 
eddy turnover rate k±5v) . We were able to reproduce 
causal critical balance of strong MHD turbulence (in the 
balanced case) by argumentation that involved field wan- 
dering. We call this effect, that comes from field wan- 
dering, the propagation critical balance and argue below 
that it will provide qualitatively new picture in the im- 
balanced case. 

Let us now consider the imbalanced case. We will ad- 
here to the same "eddy" ansatz as GS95 in that for every 
transverse scale of a wave there is a characteristic longi- 
tudinal scale. This is the physical meaning behind terms 
"eddies" or "wavepackets". We will digress from GS95 
in a natural way, postulating that both waves have dif- 
ferent anisotropies, i.e. the dependence of longitudinal 
scale A to transverse scale A is different for each kind of 
wave. This situation is presented in Fig. 1, where some 
arbitrary longitudinal scale A~ corresponds to the two 
different transverse scales, Ai for weak wave w~ and A2 
for strong wave w + . A + is a longitudinal scale of w + 
wave having transverse scale Ai. Since we originally de- 
cided to consider strong turbulence, let us assume that 
at least the w~ is being strongly cascaded by w + . In 
this case the most effective mixing of w~ on scale Ai will 
be obtained through w + motions that are on the same 
scale 2 . The longitudinal scale for w~ will be provided by 
causal critical balance, since its cascading is fast. 

2 This is the assumption that is often employed in a strong tur- 
bulence, such as Kolmogorov model of incompressible hydro tur- 
bulence, or GS95. 
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The cascading of w + is somewhat more complicated. 
Since the amplitude of w~ is not large enough to provide 
strong perturbations in w + , the w + will be perturbed 
weakly, and the cascading timescale will be diminished 
according to the "strength" of the w~ , just like it does 
in weak turbulence. Moreover, now the w~ eddies will 
be cascading w + eddies with similar longitudinal scales, 
which is the generic feature of weak cascading. 

The perturbations provided by w~ will have a trans- 
verse scale of Ai. In other words, the energy of w + will 
be transferred between A2 and Ai. What determines the 
longitudinal scale for cascaded w + ? This is the central 
question of this paper. We argue, that this longitudi- 
nal scale will be determined by the propagation critical 
balance, in the following way. The wavepackets of w~ 
are strongly aligned to the mean field on scale Ai, there- 
fore they are randomly oriented with respect to the mean 
field at a larger scale A2 . The RMS angle of wavevector 
of w~ eddies with respect to mean field on A2 will be 
around 9 « 5b + (\2)/vA- This slant of w~ wavepackets 
will determine the increase of k\\ for newly cascaded w + 
packets at Ai (see Fig. 1). 

Finally, one can verify that the increase of k\\ of w + due 
to causality is smaller (due to the fact that cascading is 
rather slow), and can be neglected. Also the increase of 
fc|| in w~ eddies due to propagation effect will be negli- 
gible due to very weak wandering provided by w~ field. 
In this picture we expect the w~ wavepackets on scale 
Ai to closely follow mean field lines on scale Ai. Essen- 
tially, this will require A~ > A + . Since the amplitude of 
the w + is much larger than that of the w~ , it provides 
much stronger increase of fc|| by propagation, effectively 
making the w + eddy less anisotropic than w~ eddy and 
validating this assumption. 

Summarizing the result of this section, the new inter- 
pretation of critical balance in the strongly imbalanced 
case is that the k\\ of the weak wave increases due to the 
finite lifetime of the wave packet, while in the strong wave 
it increases due to the field wandering of the strong wave 
itself on larger scales. This effect does not contradict 
the exact MHD solution of the wave propagating in one 
direction, because it requires the oppositely propagating 
wave as an intermediary. The de-alignment of cascaded 
strong wave is possible because the weak wave, acting as 
cascading agent, is strongly aligned with the field lines 
on scale which is different (smaller) than the scale of the 
strong wave it is acting upon. 

3. PHENOMENOLOGY 

We proceed with phenomenological treatment of imbal- 
anced turbulence using energy cascading rules for each 
wave and the new rules for anisotropies developed in pre- 
vious section. We denote Elsasser wave amplitudes in a 
sense of RMS values as w + , w~, where "+" corresponds 
to the strong wave and "— " to the weak wave. Longitu- 
dinal and transverse scales are A and A correspondingly, 
and e^ 1 are energy fluxes for both waves. 

Since we have different anisotropies for both waves, we 
are dealing with 4 quantities which are w + (X), w~(X), 
and A + (A), A~(A), so we will need 4 equations to solve 
for these functions. There are two equations for the con- 
stant flow of energy in fe-space and other two equations 
that determine frequency uncertainties, or As. 



As was explained in §2, the longitudinal scale for weak 
wave will be defined by causality, or, 



A (Ai) = v A 



w + (Ai) y 1 

Ai 



(1) 



On the other hand, the strong wave's A will be defined 
by propagation, as 



/A+(A*)y 1 _ w+(X 2 ) 



VA 



(2) 



(while §2 arguments assume A* = Ai, it could actu- 
ally be between Ai and A2, see discussion below). The 
weak wave energy flux will be determined by the strong 
shearing of the strong wave as 



Ai 



(3) 



while for the strong wave the energy flux will be de- 
creased by the wave perturbation strength as 



e+= ^ + (^-(^) . "'-(MA-(^) . /(Al/A2) . (4) 



Ai 



v A Xi 



Here /(A1/A2) is a factor which will account for the 
fact that our weak cascading is somewhat nonlocal and 
could also include the dependence on spectral slopes (see 
Galtier et al 2002), though, for the sake of simplicity, in 
this section we will assume / = 1. 

Let us now discuss a self-similar solutions of (1-4) 
which we seek in the form of 

w ± (X) = w±(X/Lr±, A±(A) = Aq (A/i)' 3 *. (5) 

The scaling exponents are calculated as (3+ — [3- — (3, 
a + = 1-/3, a_ = (1 — a+)/2, here (3 is arbitrary. The 
correspondence with GS95 is obtained with [3 = 2/3. 

Let us discuss the applicability of these solutions to 
the realistic imbalanced turbulence. Quite surprisingly, 
the self-similar solution can not be applied to the imbal- 
anced turbulence driven isotropically on the outer scale. 
This will require Aj(L) = Aq(L) which is not possible, 
because equations (3), (4) and (1) will give e + /e~ = 1. 
This situation is exacerbated by the fact, that equations 
(1-4) are nonlinear, so it is not possible to combine sev- 
eral known solutions to satisfy boundary conditions. 

Suppose, however, that self-similar solutions are real- 
izable with some outer scale boundary conditions. The 
way to uphold correspondence with GS95 (3 = 2/3 is to 
take critical balance (1) and (2) at outer scale for weak 
and strong wave respectively and modify (2) such as both 
conditions are in agreement. This is achieved, e.g., by 
replacing A + with geometrical average of A + and A - . 
Note that (3 — 2/3 is also the one that has corresponding 
spectral slopes a ± equal to each other and a Kolmogorov 
value of 1/3. If (3 < 2/3 the spectra are meeting and the 
self-similarity is broken. 

With the solution, described above, we have Aq" / A^ = 
e + /e~ and A 2 /Ai = (e+y ) 3 / 2 . Note that this solution 
is realized not necessarily due to a different anisotropy of 
w + and w~ driving, but in a quite more general setting 
where Aq (L) and Aq~(L) play the role of parameters of 
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the asymptotic power-law solution. Indeed, if we drive 
turbulence with the same anisotropy on outer scale, there 
will be a non-power-law section of the solution that will 
ensure that A 2 /Ai > 1, to fulfill e + /e~ > 1, which, given 
enough inertial range, will transit into the asymptotic 
power-law solution. We believe that this situation is re- 
alized, e.g., in the solar wind, where the counter- waves 
are generated by reflection on density inhomogcncities on 
large scales and thus have the same outer scale anisotropy 
as the direct waves. 

4. FACTORS, RELAXATION TO THE STEADY STATE, 
VISCOUS SCALE AND LIMITING CASES 

The equation (4) is written, assuming that cascading 
is weakened in a way similar to how it is weakened in the 
weak Alfvenic turbulence. This assertion does not come 
from any analytical argument, since our turbulence could 
not be treated perturbatively. Due to the fact that we 
have two length scales, Ai and A2 and two amplitudes, 
w + and w~ there are several dimensionally correct ways 
to write dissipation proportional to the fourth order of 
amplitude. We stuck with expression (4), as it seemed 
most plausible. On the second note, the cascading in 
weak turbulence is nonlocal in a sense, that the dissi- 
pation rate depends not only on the local values of per- 
turbation, but on the slope as well (Galtier et al 2002, 
Lithwick & Goldrcich 2003). Furthermore, in our model, 
the cascading of w + is somewhat nonlocal, if Ai ^ A2. 
The energy is being transferred not in the local vicinity 
of one particular scale, but between two different scales. 
In principle, this would require a logarithmic correction 
factor such as / = 1 + /o log(A2/Ai). 

Strong local driven turbulence have steady-state relax- 
ation timescales which are usually equal to the dissipa- 
tion timescale. Not only the w + cascading time is much 
larger (see (4)), but there is some nonlocality and slope 
dependence that, in principle, could carry information 
upscale and make the establishment of the steady state 
even longer. With nonlinear dissipation timescale being 
very large, the linear damping or viscosity will become 
important earlier, diminishing inertial range. These fac- 
tors make numerical study of imbalanced turbulence ex- 
tremely challenging. Also, in a realistic compressible tur- 
bulence the parametric instability or steepening become 
relatively more important, as the turbulent cascading 
rates decrease. 

Let us now discuss the onset of viscosity at small scales. 
If our turbulence would be completely local both in sense 
that same scales of "+" and "— " waves cascade each 
other and in a Kolmogorov sense that energy is trans- 
ferred locally, that would mean equal viscous scale for 
both waves and that the ratio of amplitudes must satisfy 
w + /w~ = (e + /e - ) 1 ' 2 . Our model, however, is some- 
what nonlocal, and, according to our picture of cascad- 
ing, presented on Fig. 1, strong waves w + on scale A2 
are cascaded by weak waves m" ona smaller scale Ai. 
Once the viscous dissipation becomes important for w~ 
on Ai scale, the amount of w + waves, cascaded from 
A 2 to Ai will strongly decrease. Therefore, the w + will 
have shorter inertial interval than w~ . This effect is 
observed on numerical spectra on Fig. 2. Since the 
viscous scales are different for w + and w~ , the ampli- 
tude ratio mentioned above will be offset by some power 
of Ai/A 2 . Designating viscous scale as A^ and assum- 




FlG. 2. — The power spectrum of energies for both wave species 
in an imbalanced forced incompressive 320 3 numerical simulation. 
The energy input ratio was about 4. 




Fig. 3. — The A(A) dependence in the same simulation as in 
Fig. 2, the length scales are in the units of the cube size. The A <-> A 
correspondence is established by the transverse and longitudinal 
second-order structure functions. 



ing A+ = \ v (A2/A1) 1 / 2 we obtain w+(A+)/w {\ v ) = 
(e+/e-) 1 / 2 (A 2 /A 1 ) 1 /2. if we use the self-similar, GS95- 
like solution from §3, i.e. A 2 /Ai = (e + /e~) 3 / 2 , we get 
w + (A+)/w~(\~) = (e + /£~) 1,25 (there is no guarantee, 
however, that a self-similar solution will be established 
in a limited inertial range). 

The limit of vanishing weak wave is understood as fol- 
lows. When e~ goes to zero, the steady-state relaxation 
time goes to infinity, in other words, the steady state is 
never established and the amplitude of the strong wave 
grows indefinitely from driving. Also, as A+ increases and 
approaches the outer scale, the viscous dissipation be- 
comes more important than the nonlinear dissipation and 
the cascading stops, with the strong wave being present 
only on outer scale and being dissipated viscously. 

Finally, we would like to present our preliminary nu- 
merical results of imbalanced driven turbulence. Both 
waves were driven independently with the same isotropy 
on the outer scale and different amplitudes. The energy 
input ratio e + /e~ was about 4. We evolved the simula- 
tion for more than 30 Alfvenic times unless energies of 
both species, their spectra and anisotropics has become 
stationary. The existence of the stationary state of imbal- 
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anced forced turbulence is one of the main assumptions 
in this paper as well as in LGS07. The resulting spectra 
and anisotropies are shown in the Figs. 2 and 3. The 
ratio of the total |w + | 2 to |w~| 2 "energies" was around a 
100. In the middle of the inertial interval the ratio of the 
power spectra was around 70 (LGS07 is predicting this 
ratio as (e + /e~) 2 « 16). The anisotropy curves (Fig. 
3) diverge from outer scale, as we suggested in §3. The 
w~ perturbations exhibit extreme anisotropy on small 
scales, while w + anisotropy is less pronounced. The in- 
ertial range of w~ is somewhat larger than that of w + 
(Fig. 2), which is consistent with our model above. The 
effect of pinning, i.e., the equality |w + | 2 = |u>~| 2 at the 
dissipation scale, predicted for the weak imbalanced tur- 
bulence by Lithwick & Goldreich (2003) is not observed 
in our simulations either 3 . 

At present, we think that different anisotropies of w + 
and w~ perturbations, and, most notably, the weaker 
anisotropy of w + , is the best available evidence favoring 
our model. Note that the weaker anisotropy of w + is 
hard to explain in the framework of GS95 where parallel 
scale is equal to the product of the cascading timcscale 
and the Alfven speed. 

We will present more extensive set of imbalanced tur- 
bulence simulations and a careful comparison with exist- 
ing models in a future paper. 

5. DISCUSSION 

The first treatment of imbalanced case of strong 
anisotropic turbulence, LGS07, uses subtle argumenta- 
tion in order to support the notion that anisotropies for 
both waves will be the same. The main argument was 
that even though the strain rate of the weak wave will 
be smaller, it will be applied coherently, thus resulting 
in large perturbations in the strong wave. Therefore, as 
LGS07 argued, the old critical balance will give the same 
longitudinal scale for both waves. We feel that this ar- 
gument is somewhat forced. Indeed, the coherent appli- 
cation of the strain does not necessarily mean cascading. 
While in case of the strong shearing there's only one lon- 
gitudinal scale, in the case of weak shearing there are 
two: the correlation scale and the coherence scale. This 
is clearly seen in the case of weak Alfvenic turbulence 
when the correlation scale is essentially determined by 
the outer scale and is not changed by cascading, while 

3 Our ratio, however, is beyond the applicability of Lith- 

wick & Goldreich (2003). 



the coherence scale is determined by nonlinear cascad- 
ing and is typically much larger than correlation scale 
(in fact this inequality is the condition of applicability 
of weak turbulence). In weak turbulence the large co- 
herence scale does not lead to any unusual increase in 
cascading. It is not at all obvious why this situation 
should change in the case of strong turbulence, consid- 
ered in LGS07. Another shortcoming of LGS07 model is 
that, while being completely local in both k\\ and k±, it 
dictates different nonlinear cascading timescales for "+" 
and "— " waves, which leads to inconsistency near the 
dissipation scale. Indeed, LGS07 has w + /w~ = e + /e~, 
while the local model must have w + /w~ = (e+ /e - ) 1 / 2 
to ensure transition into viscous/resistive dissipation (cf 
§4). 

We note, that earlier closure models for isotropic MHD 
turbulence used phcnomenological "relaxation of triple 
correlations" (see, e.g., Pouquct, Frisch, Leorat, 1976) 
with Alfvenic timescale, which effectively lead to the 
weakening of interaction and the IK —3/2 spectrum. 
This approach introduced the idea that large scale field 
can lead to wave decorrelation and the increase of wave 
frequency and k\\ . In this paper the increase of fey is ob- 
tained in a different way, namely, it is not dictated by 
the full background field B and the Alfvenic timescale 
(ki>A)~ j that is derived from it, but by a perturbation 
(5B, which describes a difference between local field av- 
eraged on two different scales, Ai and A2 (see Fig. 1). 

The influence of strong turbulence field wandering on 
the k vector of the wave was proposed in Farmer & Gol- 
dreich (2004) in the context of the cascading of quasi- 
parallel modes, where it was necessary to obtain non-zero 
fcj_ in order the wave to dissipate (see also Beresnyak & 
Lazarian, 2008). In our problem the situation is, in a 
sense, reversed, as we need to obtain the increase of k\\ . 
In both cases the decorrelation seems to be a feature of 
strong turbulence that does not directly appear in per- 
turbativc calculations. 
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